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Abstract - For the fractional action cosmological model, derived earlier by the author 
from the variational principle for a fractional action functional, the exact solutions 
are obtained. The case of a quasi - vacuum state of matter that fills the universe 
£\j | is considered. Moreover, on the basis of specific ansatz proposed in this paper for 

the cosmological term, the class of exact solutions of the model equations is obtained. 
" Examples for some given laws of the cosmological term evolution are provided. Besides, 

s ! , a formula for the effective equation of state is derived, and the deceleration parameter 

j^jfj of the obtained models is studied. 
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Tj'. 1 Introduction 

cr 

5^1 | As well known, the fractional action cosmology (FAC) is the new field of research based on 

the principles and formalism of the fractional calculus applied to cosmology. In this paper, 
we present the results of study of some exact models in FAC constructed by the author 
£S| ■ on the basis of the fractional action functional in pQ. These models are followed from the 

^ 1 fractional variational principle for the dynamical field theories in general and the theory of 

gravity formulated by El - Nabulsi (see, e.g., [2], [3] and references therein). In this approach, 
the action integral Sl[q] for the Lagrangian L(t, q{r), q(r)) is written as a fractional integral 



in 
rn 

O" S L [ qi \ = -^ I L(r, qi (T), qi (T))(t-Tr-'dT, (1) 



which is essentially the Riemann-Stieltjes integral for the function L at fixed t with an 

integrating function Ot(r) = — ; -\t a — (t — r) a ]. This function has the following scaling 

T(l + a) to to 

property: g fl t(l J ' T ) — ^ a 9t{ T ), fJ-> 0- It is worth to note that recently Calcagni [5], [5] gave 
a quantum gravity in a fractal universe and then investigated cosmology in that framework. 
That theory is Lorentz invariant, power-counting renormalizable and free from ultraviolet 
divergence. The action in this model is also a Lorentz-covariant and is equipped with a 
Stieltjes measure, but the model equations are somewhat different from the basic equations 
of FAC. The holographic, new agegraphic and ghost dark energy models in the framework 
of fractal cosmology are investigated in [7]. The fractal universe in which dark energy 
interacting with dark matter is considered in [8]. 

In the present work, we aim to develop further the FAC model proposed by the author 
in [T] and obtain the new exact solutions to this model. Exact solutions of the modified 
Friedmann equations in FAC are important to analyze the possible behavior of these models 
in aspect of the accelerated cosmological expansion in aspect of the recent observational 
data. To the best of our knowledge, until now, a small number of such solutions is known. 
Moreover, these solutions are mainly based on a predetermined evolution law of the scale 
factor [5] , [TU] , although it would be logical to start with some reasonable physical or math- 
ematical proposals, and only then to analyze the behavior of the scale factor. In this paper, 
we obtain several exact solutions from the assumption of a vacuum-like state of matter or 
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from a certain, rather general, relation between the cosmological term and the Hubble pa- 
rameter, which were previously used by many authors in some particular cases. Besides, 
we derive a formula for the effective equation of state (EoS) in FAC. Graphical illustrations 
of all obtained models are also given. Our models demonstrate the new types of cosmic 
evolution, which are not standard for the cosmological models, based on General Relativity. 



2 The model equations of FAC 

Following the definition (flj, the modified fractional effective action in a spatially flat 
Friedmann-Robertson- Walker metric, ds 2 = N(t) 2 dt 2 — a 2 (t)8ikdx l dx k , where N is the lapse 
function and a(t) is a scale factor, is represented by a fractional integral as follows pQ: 



a 2 dN 



o 



IttG \ N 2 N 2 N 3 




(t-T)°- l dT, (2) 



where a G (0, 1), the over-dot denotes differentiation with respect to time t, and C m is the 
Lagrangian density of matter characterized by the energy density p and pressure p. The 
variation of the action ©, followed by the choice of gauge N = 1, yields the modified 
equation of continuity, 

P+3{H+^^j (p+p)=0, (3) 

and the set of Euler-Poisson equations, which can take the form of the following modified 
Friedmann equations [T]: 

p = 3H 2 + 3 (1 ~ a) g - A , (4) 

p = -2H 3H 2 2 { -^H a) j 2 - a) +A, (5) 

t t 2 

where the gravitational constant 8irG = 1, and H(t) = a/a is the Hubble parameter. It 
is known that in standard cosmology the continuity equation for a perfect fluid, that is 
the energy conservation law for matter, follows from the Bianchi identity for the Riemann 
curvature tensor. In other words, in the standard cosmology where a = 1, the continuity 
equation ^ is a differential consequence of the field equations ((H) , (O . It is easy to prove 
that these equations also yield the modified continuity equation in the case a nel , but 
only if the following equation is valid PQ : 

g + 3^-g^ g -( 1 - Q )f- a ) = J^. (6) 
t t l 1 — a 

Note that equation © is written down after dividing by a nonzero factor (1 — a), which 
is originally appeared on the left hand side of this equation. Therefore, in the limit of the 
standard cosmology of General Relativity, i.e. for a = 1, the original equation ([5]) uniquely 
leads to a constant cosmological term A, and the set of equations (jl]), (J5J) takes its usual 
form of the Friedmann equations. 

The continuity equation (|6]) is explicitly integrable for a perfect fluid with the EoS p = wp 
with w = constant, which gives: 

P ~ 3(l + u;) t (l + u>)(l-a) ' ( ' 

If the constancy of the barotropic index w is not required as a condition, that is supported 
by the recent observational data and theoretical studies , then it is more convenient to deal 
with the set of equations ([J} - © as the main system of dynamical equations for our model. 
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Then according to the expressions dU and ([5]), the effective EoS w e ff = p e ff/p e ff, where 
p e ff = p — A and p e // = p + A, can be reduced to the following form: 



H I- a (l-a)(2-a) 
^f f = -1 - 3 1-q ' ( 8 ) 

2 ij 

It is easy to note that w e ff coincides with the known standard expression w e ff = — 1 ^ 

in the limit a — > 1, but can significantly differ from it in the case of fractional order of the 
effective action @. Moreover, the possibility of vanishing of the numerator in © at some 
instant means that the model admits crossing of the so-called phantom divide w e ff — — 1- 
Of course, the latter does not mean that the EoS of matter is necessarily w = — 1 at the 
moment of crossing the phantom divide. It should be noted also that the effective EoS is a 
dynamical characteristic of the model and gains a new definition represented by the formula 
©, but the deceleration parameter, 

a 2 a H 

is defined just as in the standard cosmology, being a kinematical parameter of the model 
[11]. 



3 The quasi- vacuum EoS of matter: w = — 1 

Let us consider a simple example of exact solution for our model which conforms to the 
following condition on the EoS of matter: w = — 1. It follows from (J7J that p(t) = po = 
constant and p = —p = —po, exactly as in the standard cosmology of GR. Then the 
rest independent equations of the system ((H) - (j6|) for the Hubble parameter H{t) and 
cosmological term A(t) can be rewritten as follows: 

Equation (fTOj) can be easily solved, and this solution for the Hubble parameter shows 
that it varies with time as follows: 

1 - a 

H=^+H t~^~ , (12) 

where C a = — — — , and H is a positive constant of integration. The latter yields 

(3 -a) 

the scale factor as a function of time t: 



3 — a ' 



a 



a = a Q & a exp ——H t 2 . (13) 



As an illustration, the graphs of functions H(t) and a(t) for a = 0.8 are shown in Fig. 1. 

The effective cosmological term A e / / = A + po as a function of time can be found from 
the equation (fTTj) in the form 

(1 + a) 

a | o ff (l-a)(7-3a) , 3(1 - a) 2 (2 - a)(5 - 2a) 1 
A e// _3ff i +3H t 2 + ^—^ -. (14) 
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-(a) B = 7.5; n = 3.5 (b) B = 8.5; n = 4.5 



Figure 1: The evolution of the effective 
equation of state w e f / , the deceleration pa- 
rameter q, the scale factor a and the Hubble 
parameter H in the model w = — 1. Here 
a = 0.8, < = 0.25 and H = 0.15. 



Figure 2: The evolution of the Hubble pa- 
rameter (j2"3"|) and the scale factor (|24[) for 
different values of i? and n. In both cases, 
A = 1, C = 3 and a = 2. 



It is seen that in the limit a — > 1 , which leads to the standard cosmology, our solution (Ti~2"j) - 
(fl4|) tends to the usual exponential expansion of the universe: a = aae^ ^ , = Hq, A e // = 

An interesting feature of the obtained solution is that according to the formula ([8]) and 
equation (JTOJ) the effective equation of state coincides with the equation of state of matter 
that is equal to w e ff = —1, but the cosmological term evolves according to the formula 
(fl"4"|) . Even more significant difference of this model from the standard ACDM model is the 
behavior of the deceleration parameter. From the equations © and (JTDJ) , it follows that 

l_a (l_aK2_a) 
q{> 2(tH) + 2{tHf ■ 

Using the explicit form of expression tH = C a +Ho £( 3 ~ Q )/ 2 from (|12p . it is easy to show that 
the deceleration parameter starts to decrease from the value go = q(t = 0) = — 1 + C^ 1 > 
down to zero, reaching it at some instant ti, which is easy to find from the equation q{ti) = 
as 

2 _ 2 

/. tt n — o r r, r (1 — — 5a)" 

ti = (4ff ) 3 ~ a V I"" 17 -9a ^ 

L 6 — a 

After that, initially slow expansion is replaced by an accelerated expansion for all t > ti. 
At time tds = [2(2 — a)/H (3 — a)] 2 /( 3 ~ a ) the deceleration parameter crosses the de Sitter 
line q = —1 and then reaches its minimum q m i n = —1 — [(1 — a)/8(2 — a)] at time t m i n = 
[(2 — a) (5 — a)/Ho(3 — a)] 2 ^ 3_Q ' > tds- From this instant, the deceleration parameter 
asymptotically tends to qoo — —1, staying in the domain q < — 1. The plots of q and w e ff 
versus time are shown in Fig. 1. 

For the standard ACDM model, i.e. when a = 1, we find ti = 0, which corresponds to 
the accelerated expansion at any time t > 0, accompanying by the constant deceleration 
parameter q = — 1. 



4 Models with a given A(t) 

Let us obtain a class of exact solutions of the equation © , assuming that the cosmological 
term A(t) is given by some condition, and it is a known function of time. Then we make the 
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Figure 3: The deceleration parameter q(t), Figure 4: The effective EoS u> e //, corre- 
corresponding to the cases (a) and (b), rep- sponding to the case (a) in Fig. 2 as a 
resented in Fig. 2. function of time for different a. 



following substitution in the equation (|6]): 

x = hx(t/t ) t = t a exp(x); Y(t) =tH(t), 
where to > is a constant. As a result, it can be rewritten as 

Y'(x) - (9 - 2a)Y(x) + 3Y 2 (x) - (1 - a) (2 - a) = -^2_e 2x A'(x), 

1 — a 



(15) 



(16) 



where the prime denotes the derivative with respect to x. From the form of this equation, we 
can assume that there exists a class of solutions for this model, for which the cosmological 
term satisfies the following equation: 



A'(x) = l —^e- 2x (k x Y\x) + k 2 Y(x) + k 3 Y 2 (x) + k 4 ) , 



f 2 
'o 



(17) 



where fcj are arbitrary constants. Indeed, substitution of p7|) in the equation (fT6"|) leads the 
latter to the form: 

AY'(x) — BY(x) + CY 2 (x) = D, (18) 



where the coefficients are: 

A = 1 - fci, B = 9 - 2a + k 2 , C = 3 - jfe 3l D = (1 - a) (2 - a) + jfe 4 . 



(19) 



Therefore, all models in which the cosmological term would be obtained from the ansatz 
(fT7|) have the same type of solution and evolve similarly. The general solution of (fT8"|) can 
be written as 



Y(x) 



(20) 



where B 2 +ACD > and C ^ 0, and the constant of integration is absorbed by the arbitrary 
constant to, which can always be done in view of the definition (|15[) . 

Note one more feature of the model, associated with a particular type of symmetry of 
the equation (fT8"|) . It is easy to see that X(x) — 1/Y(x) satisfies equation of the same type 
as (fl8|) with replacement of the coefficients p9|) according to 4 -> 4; B ^ -5; C f> f 
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Figure 5: The scale factor as a function of time for the model n = ino- Here, A = 6, B = 
8, C — 3 and ao = 1. 



provided C,D 7^ 0, that is the equation AX 7 (a:) + BX(x) + DX 2 (x) = C. Therefore, the 
solution for X(x) = Y^ 1 (x) can be represented as follows: 



-B + yjB 2 + WD tanh ^ ( x + x o) 



(21) 



. , , , ,. . , \JB 2 + 4CD y/WTICD 

where the constant xn is determined by the condition tanh ; xn = — . 

3 2A B 

After that, we can see that the effective equation of state © can be simply expressed in 

terms of the function X(x) of (|2ip as follows: 

1 2X' + (3 - a)X - (1 - a) (2 - a)X' 2 
- 1 + 3 l + (l-a)X — • (22) 



Returning to the original variables according to (|15[) . solution of the equation (|20|) for 
the Hubble parameter can be written as 



H(t) = —( K B + 2nAt & nh[n\n(t/t )}), n = ^ A , (23) 
which yields the following expression for the scale factor: 

a{t) = a t B/2C I cosh[n ki(t/t )} ) , (24) 



where ao is a constant of integration, n G R. The dependence -ff(i) and a(t) on time, 
represented by the formulas (|23p and (j2"4")) . is shown in Fig. 2. Graphs of the evolution 
of the deceleration parameter q(t) and the effective EoS w e ff(t) for the obtained solutions 
are shown in Fig. 3 and Fig. 4, respectively. It may be noted that for certain values of 
the coefficients in ((T8]) . i.e. the constants a and fcj in the formulas (|19[) . the late cosmic 
acceleration and the repeated crossing of the phantom divide w e ff = — 1 are possible. Such 
a behavior of the universe is supposed by the modern astrophysical observations [14] , [15] . 

The evolution of the model qualitatively changes in the case B 2 + ACD < 0, which is 
permissible due to the arbitrariness of the constants fcj in the definition ()19j) . Denoting in this 
case a purely imaginary coefficient n — ino, where no = \/\B 2 + ACD\ G K, we can write the 
expression for the scale factor similar to ((24j) , but with the substitution cosh[nln(t/to)] — > 
| cos[no ln(t/to)]|- This means that the evolution of the scale factor acquires the cyclic 
property on a logarithmic time variable x = ln(t/to). Graphs of the function a(t) for various 



G 



values of n are shown in Fig. 5. Note that the evolution of the scale factor of the cyclic 
nature is previously discussed in the literature in different contexts. 

Finally, let us consider some examples of with the phenomenological functions A(i), 
having an observational basis and widely discussed in the literature (see, e.g. [T2], [H])- To 
this end, we rewrite the ansatz (|17p in terms of the original variables using the definitions 
([TSl in the form 



A = (1 - a) + (fa + k 2 ) — + k 3 — + — 



(25) 



and consider the following cases. 

(1) Let fa — 0, where i — 1,..,4. Then from (|25p . it follows that A = constant. In 
accordance with the definitions (|T9]) and (|23| , we obtain: 



f 



A t = l, Si =9- 2a, Ci = 3, m = - V(9 - 4a) 2 + 24. 

(2) Now let fci = fc 2 = fc 3 = and fcj = —2/3/(1 — a), where /3 is an arbitrary positive 

constant. Integrating (|25|) . we obtain A(i) = where the integration constant is assumed 

t z 

we find the following coefficients: 
1 



to be zero, and from (|19l) and 

A 2 = l, B 2 =9-2a, C 2 =3, n 2 



,. / (9-4a) 2 + 24 1 a 13 . 
2 V 1 — a 



(3) Finally, let fa 



, k 2 = —2fa, k 3 = &4 = 0. From equation (|23|) it follows 



that A(t) = —Hit), and the constant coefficients of (fT9)) and (|23|) are given by the following 
expressions: 

A 3 = 1 ~ a ~ , B 3 = 9~2a 



n 3 



1-a-P 
1 - a 
1 



2(1 - a-/?) 



l-a 



, C 3 = 3, 



(1 - a) (9 - 4a) - 2/3 + 8(1 - a) 3(1 - a) - a/3 



It is interesting to note that the parameter n 2 can become imaginary for the coupling 
constant /3 > (1 — a) [(9 — 4a) 2 + 24]/24, which corresponds to the case illustrated in Fig. 5. 



5 Conclusion 

Thus, in this paper we have obtained some exact solutions of the cosmological model derived 
from the variational principle for fractional action functional in [I]. Exact solutions for the 
dynamical equations of this model are found either under the assumption of a vacuum-like 
state of matter that fills the universe, or on the basis of a rather general ansatz for the 
dynamical cosmological term. In any case, the behavior of the model demonstrates a signifi- 
cant difference from the corresponding standard model, which obviously is a consequence of 
the fractional order of action functional and, possibly, the fractal nature of space-time (see 
e.g. [5], i). 

Particularly attractive in our opinion is that the effective EoS which is defined by pj). 
can cross the phantom divide line, as supported by the recent astrophysical observations. 
Recall that in standard cosmology with the sole source this crossing is not possible. In 
addition, the model is able to evolve cyclicly, that is attracting the attention of researchers 
for the recent years (see [16] for a review). 

Finally, we have to note that our ansatz (fl~7j) for the cosmological term proposed in the 
present work, can not cover all the features of the models in FAC. Therefore, it would be 
interesting to continue this research towards a more general ansatz, which could include the 
greater number of phenomenological expressions for A(<). 
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